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1. Introduction 

Models for discrete polymerized membranes with rigid bonds display interesting phys- 
ical behaviors, as the only way to change their spatial configuration is through folding. We 
will consider here only "phantom" folding, in which the membrane is allowed to interpen- 
etrate itself, as opposed to self-avoiding folding, which introduces considerable mathemat- 
ical difficulties, as was already noticed in the case of compact (one-dimensional) polymer 
folding 0. 

The entropy of folding of a regular triangular lattice has been first studied numerically 
in and later computed exactly in 0], by mapping the model onto the edge tri-coloring 
problem of the triangular lattice. The next step was the introduction of a bending rigid- 
ity, namely an energy which favors either folded edge states or unfolded ones (typically, 
one associates the Boltzmann weights e ±K respectively to each unfolded or folded edge). 
Numerical studies have shown the existence of a first order folding transition between 
completely folded states (with a certain entropy) and a unique completely flat state, at 
a critical value K = K c > 0. On the other hand, the model also displays a continuous 
unfolding transition between a completely folded state and a partially unfolded one at 
another critical value K = < 0. Only little numerical evidence for the latter transi- 
tion was found in Q, but it was argued that for K < 0, starting from a completely folded 
state, edges could be unfolded along loops, a phenomenon typical of continuous (Ising-like) 
transitions. 

In the present paper, we carry out an analogous study for the square-diagonal lattice 
folding problem introduced in ||. The square-diagonal lattice (see Fig.l below) has two 
types of edges, long and short, it is therefore natural to attach to them two different types 
of Boltzmann weights, e ±Ki , i = 1, 2 for respectively short and long edges. In addition to 
getting better numerical estimates when = or K\ = 0, we will be able to investigate 
the full critical structure of the model in the plane (Ki, K%). 

The paper is organized as follows. In Sect. 2, we review the 2-dimensional folding 
problem of the square-diagonal lattice, and its reformulation as an edge tangent-vector 
model, and various face/vertex models. In the presence of bending rigidity, the transfer 
matrix of the model is simple enough to permit numerical estimates of the free energy and 
magnetization, by extraction of its largest eigenvalue. Sect. 3 is devoted to the numerical 
study of the rigid short edge model, with K2 = 0, in which only the short edges have a 
bending rigidity, whereas Sect. 4 deals with the rigid long edge model, with K\ = 0. Both 
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models confirm the first order folding transition of [Q, at some critical values K = K c ^ > 0. 
In Sect. 5, the two above models are studied for Ki < 0, and better evidence is found for a 
continuous unfolding transition at some critical values Ki = K^* < 0. The general rigid 
long and short edge model is studied in Sect. 6, where we find the complete phase diagram 
of the model, including two critical curves corresponding to the folding and unfolding 
transitions, and which are separated by the line K<i = —2K±. 

2. Folding the Square-Diagonal Lattice 




Fig. 1: The Square-Diagonal lattice. It has two types of vertices, respectively 
4- and 8-valent, and two types of edges, short (length 1) and long (length \/2)- 

Let us consider the Square-Diagonal lattice, obtained from the standard square lattice 
by drawing the two diagonals on each face, as depicted in Fig.l. Each face of this lattice 
is a triangle with two "short" edges of length 1 and one "long" edge of length y/2. In 
this paper we study the possible foldings of the lattice, in which the edges serve as hinges 
between adjacent faces. The simplest such foldings are those with a two-dimensional final 
state, namely the foldings of the lattice into itself. Here we consider only phantom folding, 
in which the lattice may interpenetrate itself, and distinguish only between distinct final 
two-dimensional (folded) states of the lattice, namely its folding configurations. 

In |3| , we have shown that this two-dimensional folding problem is equivalent to various 
vertex models. We first defined the tangent vectors to the lattice as the vectors t drawn 
along its edges, oriented in a compatible way, so that the sum 

face 
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taken around each face of the lattice vanishes. The short tangent vectors may take the 4 
values ±e~i, ±e~2, where (€{,€2) is an orthonormal basis for the short edges, whereas the 
long edges may take the 4 values ±e~l ± e~2. 

A folding configuration of the lattice is entirely determined by the images of these 
tangent vectors under a continuous, length-preserving folding map p, with the constraint 
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(2.2) 



face 



around each face of the lattice. Note that the images of the tangent vectors under p may 
take only values in the corresponding abovementioned sets of 4. 

By focusing on either short or long edges (the images of either class are sufficient 
to reconstruct all other images, thanks to the rule ( |2.2| )), we may consider the folding 
configurations as either short or long-edge configurations on the lattice, with constraints 
inherited from the rule ( |2.2| ). In both cases, we end up with statistical models with 4 edge 
variables, and only finitely many allowed vertices. 
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Fig. 2: The 28 possible face configurations of the rigid short edge model. 
The variable a may take any value mod 4. We have indicated the Boltzmann 
weights including a bending energy K\ per short edge and a magnetic field h. 



On long edges, the constraint (|2.2| ) amounts to relating the four long edge images 
around each square face made of 4 triangles. Let us first introduce the orthogonal basis 
$2) = (el + e~2, e~2 — e{) for long edge vectors, and label the edge values /1, /2, —fi, —fi 
respectively by 0, 1,2,3. In terms of these values, the only edge configurations around a 
square face which are allowed by (|2.2| ) are the 28 cases depicted in Fig. 2. They correspond 
respectively to the following folding states of the inner short edges: all four folded (4), 
two folded along the same diagonal (16), none folded (8). In Fig. 2, we have represented 
in dashed lines the unfolded short edges, and in solid lines the folded ones. Introducing 
a bending energy J\ to distinguish between folded and unfolded short edges, we are led 
to weight the first case of Fig. 2 by e -4 ^ 1 , the next four by 1 and the two last ones by 
e 4Kl , where K\ = —J\/kT. To complete the definition of the model, we need to introduce 
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a suitable order parameter for the folding transition. This parameter should distinguish 
between flat and folded phases. Following Q, we introduce a fictitious magnetic field H, 
coupled to the normal vectors to the faces, defined in such a way that they all point up in a 
completely flat configuration of the lattice. The corresponding magnetization is the order 
parameter for the folding transition. Setting h = —H/kT, the two last face configurations 
of Fig. 2 receive extra Boltzmann weights e ±4,t respectively, as one of these two flat faces 
points up and the other points down. This model is studied in Sect. 3 below. 
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Fig. 3: The 32 possible face configurations of the rigid long edge model. 
The variables a, b may take any values mod 4, with a ^ b mod 2. This gives 
8 configurations for each of the first cases, and 4 for each of the second and 
third ones. We have indicated the Boltzmann weights including a bending 
energy K2 and a magnetic field h. 

On short edges, the constraint (|2.2| ) implies first that the two short edges of each 
triangle must have perpendicular images. Using again a Z4 short edge variable reading 
respectively 0, 1, 2, 3 for the edge values el, e~2, —el, — e~2, this means that a 7^ b mod 2 for 
the two short edge values a, b of a given triangle. The constraint also relates the short 
edge images of any two triangles sharing a long edge as follows. If u, v stand for the 
short edge images of the first triangle, and u' , v' the second, there are only two ways to 
accommodate u + v = u' + v', namely u' = u and v' = v, or u' = v and v' = u. The 
corresponding configurations are depicted in Fig. 3, for the two types of faces (according 
to the position of the long edge). In the first cases, the common long edge is folded (solid 
line), whereas in the second and third it is not (dashed line). Introducing a bending energy 
J2 to distinguish between folded and unfolded long edges, we decide to weight the first 
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cases with a Boltzmann weight e 2 , and the second and third ones with a weight e 2 , 
with K 2 = —J2/kT. As before, we introduce a magnetic field H coupled to the normal 
vectors to the membrane. This results in an extra Boltzmann weight e ±2h for respectively 
the second and third cases of Fig. 3, with a flat face respectively pointing up and down. 
This model is studied in Sect. 4 below. 

The general model including bending rigidity for both types of edges is studied in 
Sect. 6. In all cases, we will construct a row-to- row transfer matrix T for rows of width 
L = 1,2, ... and extract its largest eigenvalue, which dominates the partition function 

Z = Tr(T M ) ~ \™ ax „ e LM f (2.3) 

when L and M become large. This will lead to numerical estimates for the thermodynamic 
free energy / per site (note the change of sign and the absence of the usual 1/kT in the 
definition of /, for notational simplicity). All the transfer matrices we are dealing with 
are sparse due to the relatively short number of allowed face configurations, compared to 
the size of the matrices. We therefore use an efficient algorithm which codes the matrices 
as the list of their non-vanishing elements, together with their position, and then extracts 
the largest eigenvalue by repeated iteration on a given vector. 



3. Rigid Short Edge Model 

Including the abovementioned magnetic field h = —H/kT, the Boltzmann weights for 
the face configurations of Fig. 2 read 

Wl = e~ 4Kl w 2 =w 3 = w A =w 5 = l W6 = e 4Kl+4h , w 7 = e 4Kl - 4h (3.1) 

Note that the effect of the magnetic field is to distinguish the configurations in which the 
normal vectors all point up (energy —AH) or all point down (energy AH), from all other 
cases, where there are as many up and down-pointing normal vectors. 

For definiteness, let us denote by w(a,b,c,d) the Boltzmann weight attached to a 
cyclic configuration (a, b, c, d) of edge variables, clockwise around the square face with a 
in the west position. The row-to-row transfer matrix of the model for a row of L square 
faces then reads 

L 

T au ..., aL - a ' 1 ,...,a' L = Yl Y[w(b i ,a' l ,b l+1 ,a i ) (3.2) 

6i,...,& L +iez4 *=i 
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Fig. 4: Magnetization versus magnetic field for the rigid short edge model at 
K\ = 0, for strips of width L = 1, 2, 6. We see clearly two phases: folded, 
with M — > 0, and flat, with M = 1. We expect the transition between the 
two to be discontinuous in the thermodynamic limit L — > oo. 

This matrix has size 4 L x 4 L and is sparse, as there are only ~ 7 L non-vanishing elements 
(given the west value, there are exactly 7 allowed face configurations). The numerical result 
of H for the entropy of folding / = ssd — -230... is readily recovered for K\ = h = 0. 

We display in Fig. 4 the magnetization M = dhf for various sizes L = 1,2, ...,6 at 
K\ = 0. Even for these finite values of L, we see a drastic jump in the magnetization at a 
value h = h^L{K\ = 0) of the magnetic field. Extrapolating these, we find that they tend 
to a limiting value h c (0) = .23... for large L. In the \ow-h phase, the magnetization has a 
small finite slope, which tends to for large L. In the large-/i phase, the magnetization is 
identically 1. We conclude that the system is reduced to two phases, a completely folded 
one for h < h c (0) and a completely flat one for h > h c (0). Note that h c (0) ~ ssd- 

This two- phase situation is observed for other values of K\ in Fig. 5, where we have 
represented the corresponding magnetizations as functions of h for L = 6. Like in 0, this 
leads us to formulate the hypothesis that, at least for small enough K\, there will always 
exist only two phases, folded and flat, with M = and 1 respectively. Let fo(Ki) denote 
the free energy per triangle of the model in zero magnetic field. This is exactly the energy 
of the folded phase of our model, as the latter is insensitive to the magnetic field. Note 
that this phase is entropic, at least when K\ = 0, where /o(0) = ssd = .230... On the 
other hand, the completely flat phase consists of configurations made only of faces of the 
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Fig. 5: The magnetization as a function of the magnetic field for strips 
of width L = 6, and various values of K\ = .5, .25, .1, 0, — .1. A change 
of behavior is observed for K\ = .5: there, we expect the folded phase to 
disappear in the thermodynamic limit. 

form (a, o + l, a + 2, a — 1), with weight W4 = e 4 ( Kl+h \ hence resulting in a thermodynamic 
energy fi(K± : h) = K\ + h per triangle. The critical field h c (K\) is then defined by the 
equality of the two energies, namely 

h c (K x ) = /„(#!) -Ifi (3.3) 

Let Xmlx(Ki) denote the largest (Perron- Frobenius) eigenvalue of the transfer matrix 
TofQ. Taking (j> Q , L {K x ) = Log(A^L (Ki))/{AL) as a sequence of approximations 
for /o(-ftTi), the folded phase energy, is not very accurate, as it includes a generic two-fold 
degeneracy due to the invariance under reversal of all non-horizontal tangent vectors, i.e. 
a folding along the top or bottom horizontal line of the row (in which, in particular, the 
two face configurations (a + 1, a + 2, a — 1, a) and (a — l,a + 2,a + l,a) are exchanged). 
We should rather use 

fo, L (Ki) = </> ,L(^i)-^Log2 (3.4) 

by explicitly dividing out the degeneracy. Although in principle equivalent to ^o,l(-^i) 
in the thermodynamic limit, /o,l displays a much better convergence. Let us now turn 
to the finite L flat phase energy /1 £,(l£i, h). The flat phase matrix elements of T are 
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Fig. 6: The critical field h c ^{K) for L=l,2,...,7 (from top to bottom). We 
have represented in dashed line the extrapolated limiting curve for L — > oo. 
It displays a transition at K = K c \ ~ .27. 

non-diagonal, and read T a ^ +2} a,a+2,...-a+2,a,a+2,a,... = 2e 4LKl cosh(4L/i). They give the 
energy 

h, L {K u h) = K 1 + h+ ^Log(l + e~ 8Lh ) (3.5) 

Due to the exponential character of the correction, the approximation to the critical field, 
h c ,L(Ki) = fo,L(Ki) — K\ is accurate, at least in the large h phase. We have represented 
in Fig. 6 the critical field h c ^(K) as a function of K, for various values of L. These curves 
are extrapolated to the large L limiting curve represented in dashed line. This curve has 
also been obtained by the extrapolation of another sequence of approximations, namely 

h c , L/L+1 (K) = lLo g (^^l)-K (3.6) 

for L = l,2,...,6, with an excellent agreement. We observe the existence of a finite value 
K c i = .27 of Ki beyond which h c vanishes identically. 

4. Rigid Long Edge Model 

Including the magnetic field h = —H/kT, the Boltzmann weights of the configurations 
of Fig. 3 read respectively 

Wl =e- K * w 2 =e K2+2h W3 = e K *- 2h (4.1) 
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As before, let w(a,b,c,d) stand for the Boltzmann weight of the configuration (a,b,c,d) 
around the face, clockwise with a in the western position. By a slight abuse of notation, we 
will use the same letters to denote the Boltzmann weights, transfer matrix, energy, etc... 
of the present model as those used for the Rigid Short Edge model of previous section. 
With these definitions, the transfer matrix for a row of width L reads 

L 

T ai ,...,a L ;a' 1 ,...,a' L (h,K) = ^ wfc, c., b i+1 , di)w(di, a • , d i+1 , a ) (4.2) 

Indeed, because of the alter nance of the long edge between first and second diagonal in the 
short square faces, we must consider two rows of square faces, with a total of 4L triangles. 
Let us note however that T = U{h, K)U t (—h, K), where U{h,K) is the transfer matrix 
for the first row, with entries 

L 

U ai ,...,a L ;a' 1 ,...,a' L {h,K) = Y[w(b it a' it bi+ lt a*) (4.3) 

Indeed, upon transposition, the picture for U is reflected wrt a horizontal line, hence the 
inner long edge on each face switches to the other diagonal, which amounts to a change 
h — > —h according to Fig. 3. 
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1.22891 
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6.040895 
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1.24610 
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14.36690 


1.24926 
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1.25351 
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1.25482 
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86.97476 


1.25558 
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216.4118 





Table I: Numerical results for the maximum eigenvalue of the square root of 
the transfer matrix of the rigid long edge model with free boundary conditions. 
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We have represented the size n of the row, the largest eigenvalue and the ratio 
v n = (An+i/A^) 1 / 2 , which converges to the partition function per triangle. 

This remark makes the numerical study much simpler, as the extraction of the largest 

eigenvalue of T is done by iterating the action of U and U l in alternance on a given vector. 

Again, U is sparse, with a number of non-zero elements ~ 2 L much smaller that its size 

4 L . As an application, we list in Table I the largest eigenvalue of the matrix \/T = \ / UU t 

for K 2 = h = 0, and for L = 1,2, ...,10. This is in agreement with the estimate for the 

partition function per triangle zsd = 1.258... of ||. 




Fig. 7: Magnetization versus magnetic field for strips of width L = 7 in the 
rigid long edge model, for various values of K. We note the existence of two 
phases with respectively M = (folded) and M — 1 (flat). 

The two-phase hypothesis of the previous section is confirmed here in Fig. 7, where 
we have represented the magnetization as a function of h for L = 7 and various values 
of Ki- An analogous reasoning leads to the determination of the critical magnetic field 
h c (K 2 ) = fo(K 2 ) — obtained by equating the energies of the folded and flat phases. 
Here the flat phase is made of only the face configuration (a, a + l,a,a + 1), hence an 
energy of h(K 2 , h) = (2L(K 2 + 2h))/{AL) = ^ + h. 

For strips of finite size L however, the folded phase energy reads 

foAK 2 ) = ±rLog\£l(K 2 )-±rLo g 4 (4.4) 
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Fig. 8: The critical fields h c ^{K) for L = 1,2, ...,9. We have represented 
the extrapolated large L limit in thick dashed line. It displays a transition at 
the critical value K c2 — -54. 

where \max{K 2 ) is the largest (Perron-Frobenius) eigenvalue of the transfer matrix T ( |4.2j ) 
for a row of 2 x L squares. Note the subtraction of the 4 = 2x2 times degeneracy under 
the previously mentioned reversal of non-horizontal tangent vectors, in which in particular 
(a + 1, a, a + 1, a) <-> (a — 1, a, a — 1, a). On the other hand, the flat phase energy corresponds 
to the diagonal matrix elements T 0>aia) ... ;0>a)a ..., 06/4, with only face configurations of the 
type (6, a, 6, a) and b 7^ a mod 2. This gives four identical diagonal blocks, with eigenvalue 

\%l(K 2 ,h) = 2e 2LK * cosh(4Lh) (4.5) 

and the energy 

h, L (K 2 ,h) = + h + ±Log(l + e~ 8Lh ) (4.6) 

Neglecting the exponential correction, the equality between ( [4.6] ) and (|4.4| ) yields the finite 
size critical field 

h c , L (K2) = h,L{K*)-?± (4.7) 

which we have represented for various values of L in Fig. 8. We observe an excellent 
convergence of the quantity h c ^(K 2 ), and a particular value K c>2 — .54 of K 2 beyond 
which the critical field h c vanishes identically, hence only the flat phase survives. Note 
that we get K Cj2 ~ 2K Cj i. 



11 



5. Continuous Transitions for K < 

As observed in 0] in the case of the triangular lattice folding, both short and long rigid 
edge models undergo continuous phase transitions for some negative critical values 
and K* 2 of K. In terms of our initial setting, these are anti- ferromagnetic transitions. 
However, they take place in the entropic (folded) phase of both models, in which local 
excitations can be formed, as explained in ||. 

a a 



a 



a 



a a 

Fig. 9: A local excitation of one of the anti-ferromagnetic (completely folded) 
groundstates of the rigid short edge model. The unfolded short edges are 
represented in dashed line. All others remain folded. 

More precisely, in the rigid short edge model, starting from one of the (completely 
folded) anti-ferromagnetic groundstates, in which all the short edges are completely folded, 
hence where all long edges take the same value a, we can unfold a square of 8 short edges, 
as shown in Fig. 9. This is a local excitation of the anti-ferromagnetic groundstate, which 
receives a relative Boltzmann weight 2e~ 16Kl , including the effect of unfolded edges and 
the free choice of the two orientations b = a ± 1 of the new tangent vectors. Note that 
we have not specified the folding state of the long edges, which results in a non-vanishing 
entropy of the model. Such local excitations can be combined so as to form any loops of 
unfolded short edges (with the constraint that these edges can only be unfolded by parallel 
pairs within each square of long edges) . 

In the rigid long edge model, similar excitations of the (completely folded) anti- 
ferromagnetic groundstate may occur. They correspond to unfolding four long edges form- 
ing a square, as shown in Fig. 10. This excitation receives a relative Boltzmann weight 
2e~ 8K \ A gain, these excitations may be combined into any loop of unfolded long edges. 
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Fig. 10: A local excitation of one of the anti-ferromagnetic (completely 
folded) groundstates of the rigid long edge model (with a ^ b mod 2). The 
four dashed long edges have been unfolded in the process, while all others 
remain folded. 




Fig. 11: Specific heat vs K for the rigid short (1) and long (2) edge models. 
The pairs L/L + 1 used to evaluate the curves are indicated. Both models 
display a continuous transition. The critical values of K are estimated to be 
K* tl = -.08 and K^ 2 = -.1- 

In both cases, the resulting loop model is naturally expected to undergo a phase tran- 
sition of the Ising type, when Ki, K 2 < cross some critical values. To check this, we have 
represented in Fig. 11 for L = 1,2, ...,6 or 7 the specific heats C v \ } l{Ki) = d"^- (pi^(Ki) 
and C V 2,l{K2) = d] < - 2 (p2,L(K 2 ), where the approximations to the free energy are taken to 
be 

^ ■ MIS) (5A) 
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We observe the typical growth of the maximum of the specific heat with L, which 
points to critical values of K reading respectively K*^ = —.08 and if*, 2 = — •!• Using 
the graphs of Fig.ll, we have also obtained estimates for the thermal critical exponents ctj 
such that in the thermodynamic limit C v i(Ki) ~ \Ki — K*^\~ ai , with the results a\ ~ .2 
and a 2 — .25. 

6. General Phase Diagram 

Let us now address the general rigid (long and short) edge model for the square- 
diagonal lattice folding problem. This is simply a combination of the two models of Sects. 3 
and 4, where the long and short edges receive the Boltzmann weights e ±Kl (short) and 
e ±I<2 (long) for unfolded (+) and folded (— ) states of the edge. 

In view of the results of Sects. 3, 4 and 5, we expect a very interesting phase diagram 
for this general model, in which both first order and continuous transition can be found. 




Fig. 12: The row-to-row transfer matrix T for the general rigid long and 
short edge model. T acts by transferring a zig-zag row {(ai, &i), (a2, 62), • ••> (a>L, 
of pairs of adjacent short edges to the nearest vertical translate. 

The transfer matrix T of the general model reads as follows. We consider a row of 
4L triangles, as shown in Fig. 12, which transfers a zig-zag line of 2L short edges to the 
nearest translate above it. The short edges in a line naturally go by pairs of adjacent edges 
sharing the same triangle, with a total of 8 possible images, a, b G Z4, with a 7^ b mod 2, 
expressing the fact that the corresponding tangent vectors are orthogonal. Hence we may 
view T as a 8 L x 8 L matrix acting in the space of pairs of adjacent short edges along a 
zig-zag line of 2L edges. 

The matrix T is made of L cells, expressing the various Boltzmann weights and con- 
straints on the short and long inner edges. We distinguish three pieces, according to Fig. 13, 
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w i 

Fig. 13: The transfer matrix of the rigid long and short edge model is 
made of three elementary Boltzmann weights iwi, w 2 , w% corresponding to the 
arrangements of inner short edges c,d,e G Z 4 , with c ^ d mod 2, and outer 
short edges a, 6, a', 6', a" G Z 4 , with a ^ & mod 2, a' ^ 6' mod 2. 

with the Boltzmann weights 

wi{a,b,c,d) = e K2 5 a , c 5 bjd + e~ K2 5 a:d 5 bjC 
w 2 (a',b',c,d) = (e 2K ^ a ^ c+2 + e- 2K ^ a ^ c )(e 2K ^ ¥ , d+2 + e- 2K ^ ¥td ) (6.1) 
w 3 (b',a",e,c) = e K2 5 b > je 5 a » )C + e~ K2 5 a » )b >5 e , c 

where 5 Pyq = 1 if p = q mod 4, otherwise. The transfer matrix elements then read 

T {ai,bi}iW v %} = 
L 

^ J\w 1 (ai,bi,Ci,di)w 2 (a' i ,b' i ,Ci,di)w3(b' i ,a i+1 ,di +1 ,Ci) ( 6-2 ) 

i=l,2...£ 



Inspecting (|6.1| ), we see that T has ~ 16 L non-vanishing entries, for a size 8 L x 8 L , 
hence is a sparse matrix. But the size is of course much more limiting than in the cases of 
Sects. 3 and 4. 

The transition curves we are looking for are characterized numerically, for finite size 
L of the transfer matrix, by a growing peak of the specific heat. To identify those in the 
Ki,K 2 plane, we have located the maxima of the specific heat along lines through the 
origin. More precisely, taking K\ = K cos(9) and K 2 = Ksm(0), we have taken a discrete 
set of angles 9 G [—n/2, tv/2), and represented the maxima of the specific heat 

C VtL (K) = d^Log(\£l(Kcos(e),Ksm(9)f) (6.3) 
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Fig. 14: Maxima of the specific heat C Vj l(Ki, K^) for L = 2 (stars), L = 3 
(circles) and L = 4 (diamonds), in the (Ifi,!^) plane. We have represented 
the axes as well as the line Ki = — 2K\ in dashed lines. The latter appears 
to be the separation between the continuous (anti-ferromagnetic) transition 
curve and the first order (ferromagnetic) one. 

where Xml x (Ki, K^) denotes the maximum (Perron- Frobenius) eigenvalue of the transfer 
matrix T defined above, for a strip of width L. The results for L = 2, 3, 4 and 9 = /c7r/(2m), 
k = —to, — m + 1, to — 2, to — 1, to = 35, are displayed in Fig. 14. 

Note that the separation between the two transition curves, continuous anti- 
ferromagnetic (negative -fT's) and first order ferromagnetic (positive -fT's) appears to lie 
precisely on the line K<i = —2K\. This should not be too much of a surprise, if we note 
that there are exactly twice as many short edges as long ones in the square-diagonal lat- 
tice. The anti-ferromagnetic and ferromagnetic bending energies are precisely balanced 
along the line K2 = —2K\. As a check of the crossing between the two curves, we have 
represented the specific heat C V ^(K\, —2K{) for L = 4 in Fig. 15. We see clearly the four 
critical points, two of which lie on the continuous transition curve, and two on the first 
order one. 

We are naturally led to conjecture the thermodynamic picture of Fig. 16 when L — ► 00, 
with two transition curves, continuous (anti-ferromagnetic) in dashed line, and first order 
(ferromagnetic) in solid line. Outside of the solid curve, and above the line K2 = —2K±, 
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Fig. 15: The specific heat C v ^(Ki : —2K\) for L = 4, along the crossing line 
K<i = —2K\. We have indicated the transition points (or remnants thereof). 

the model is always in the completely flat phase with magnetization M = 1 if h — > + and 
M = — 1 if h — > 0~. Outside of the dashed curve, and below the line K 2 = —2K±, the 
model is in a completely folded phase, with possibly some unfolded loops of edges. The 
staggered (anti-ferromagnetic) magnetization M s t corresponding to a staggered magnetic 
field h s t (taking alternating values ±h s t on adjacent triangles) is a function of the K J s 
tending to 1 if h s t ~~ > + and —1 if h s t — > 0~ for large negative K's. The dashed curve 
corresponds to the vanishing of the staggered magnetization, which becomes identically 
zero inside the curve. So the region inside both curves corresponds to both magnetizations 
equal to zero: we can think of it as a partially folded (or unfolded) phase. 

To fully illustrate the ferromagnetic transition, we have represented schematically 
(most of the curves are represented as lines) in Fig. 17 the surface h = h c (Ki, K 2 ) of the 
critical magnetic field beyond which the membrane is in the completely flat phase, and 
below which it is folded. This surface is singular along the plane K 2 = —2K\\ h{K\,K<i) 
has a jump discontinuity from one of its sides to the other, except at K\ = K% = 0. The 
corresponding point h = ssd = h c (0, 0) is a conic singularity of the surface. 
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K 2 — 2Kj 



Fig. 16: The extrapolated critical lines of the short and long rigid edge 
model in the plane {KijK?). The solid line curve represents the first order 
(ferromagnetic) critical line, involving in particular values K\ > and K<i > 
0. Large positive Ki's and -ftVs correspond to ferromagnetic order, i.e. the 
completely flat phase. The dashed line curve represents the continuous (anti- 
ferromagnetic) Ising-like transition line, involving in particular values K\ < 
and K2 < 0. Large negative K\& and correspond to anti- ferromagnetic 
order, hence a completely folded phase. 
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Fig. 17: A rough sketch of the surface h = h c (Ki, K 2 ) for the critical 
magnetic field separating the flat and folded phases of the membrane. This 
surface coincides with the plane h = outside of the transition curve and 
in the domain K 2 > —2K\ Note the conic singularity at the point (0, 0, h = 
ssd), and the jump singularities along the plane K 2 = —1K\. 
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